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Abstract 


Potentially serious environmental effects wxlst when cargo orbital 
i ansfer vehicle (COTV) ion propulsion is used on the scale proposed in 
the preliminary definition studies of the Satellite Power System. These 
effects of the large scale injections of ion propulsion exhaust in the 
plasmasphere and in the outer magnetosphere are shown to be highly model 
dependent with major differences existing in the predicted effects of two 
models--the ion cloud model and the Ion sheath model. The expected total 
number density deposition of the propellant Ar + in the plasmasphere, the 
energy spectra of the deposited Ar + and time-dependent behavior of the Ar + 
injected into the plasmasphere by a fleet of COTV vehicles differ drastically 
between the two models. The major environmental effect of the former model 
is communication disturbance due to plasma density irregularities (Curtis 
and Grebowsky, 1980b) , in contrast to the spectacular predictions of the 
latter model which include power line tripping and pipe line corrosion 
(Chiu et al., 1979b). The ion sheath model is demonstrated to be applicable 
to the proposed Ar + beam physics if the beam is divergent and turbulent 
whereas the ion cloud model is not a realistic approximation for such a 
beam because the "frozen-field" assumption on which it is based is not valid. 
Finally, it is shown that the environmental effects of ion propulsion may 
be mitigated by the appropriate adjustment of the beam parameters. 


Introduction 


In this initial assessment we have concentrated cn understanding 
the physics of energetic ion beam propagation in the terrestrial magnet < 
sphere. Such an understanding is essential in order that the potent i ll 
environmental effects of the projected massive amounts of ion beam 
injections during Satellite Power System (SPS) construction may be 
evaluated. Although experimental ion beam injections in space are 
required for a full understanding of the beam physics, we believe that a 
realistic beam propagation model consistent with the underlying physics 
is needed. Toward this end we will compare and contrast two proposed 
models for the physics of the ion exhaust ejected by the cargo orbital 
transfer vehicle (COTV), namely: our ion beam sheath loss model (Curt 's 
and Grebowsky, 1980a) and an ion cloud model used originally to describe 
barium release experiments in the magnetosphere (Scholer, 1970). 

The possible environmental effects fr the two proposed beam 
propagation models are radically different. For the ion beam sheath 
model we would expect a depletion of the high energy ions in the radiation 
belt (Curtis and Grebowsky, 1980b) whereas for the ion cloud model the 
trapped relativistic electrons would increase (Chiu fit al., 1979b). The 
expected radiation belt dosages thus go in opposite directions depend in; 
on the model chosen. In addition, the extension of the radiation belt .s 
out to the geosynchronous orbit with the associated radiation dosage 
hazard for the Satellite Power System (Cladis and Davidson, 1980) is 
only possible in the ion cloud model with the much more massive deposi- 
tion of thermal exhaust ions it predicts. Also, the generation of 
artificial ionospheric currents is a direct product of only the ion 
cloud model. If the assumption of frozen-in field lines does not hold 
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(which is to be expected in e manifestly turbulent pleeae environment 
■uch e* the BxheuBt ion been) , the velidity of the ion cloud model and 
ita spectacular envlromentel effecta auch as power-line tripping and 
pipeline corrosion (Chiu at al. , 1979b) would vanish. The major environ- 
mental effect that is expected from the ion beam shaath loss model that 
we have developed is the degradation of satallita communications; ground 
to space , space to ground and spaca to spaca. These communication dis- 
ruptions are caused by small spatial irregularities produced by: (1) 
plasma Instabilities in the plasmasphera driven by the energetic aniso- 
tropic ions shad by the beam during its propagation through the plasma- 
sphere; and (2) by the precipitation of magnetospheric ions into the 
ionosphere by the shed beam ions. 

Due to the drastic differences in the environmental imparts predicted 
by the two beam models, it is imperative that it be determined which of 
the models is consistent with the physics of the proposed COTV lou beam 
configuration. We will show that only the ion sheath loss model can 
adequately explain the physics of the beam. Later, we also discuss 
possible experiments to resolve beam physics and strategies that can be 
used to reduce environmental effects which may result from the beam 
injections. 
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Model Comparisons and their Related Env lronmental Effects 

In lifting meislve space power system peyloade from low earth orbit 
to geosynchronous earth orbit. Cargo Orbital Transfer Vehicles (COTV) 
using ion propulsion will Inject energetic beams of argon ions Into the 
plasmasphere (Department of Energy, 1978). The characteristic speeds 
of the plasmasphere relevant for ion beam propagation are the Alfven 
speed and the Ion thermal speed. Near the earth's equatorial plane, the 
argon ion beaus have a speed V^ < V^, the plasmasphere Alfven speed, and 
V b ^ V th* the plasmasphere ion's thermal speed (see Figure 1). The 
Alfven speeds are shown for the Chiu at al. (1979a) model plasmasphere 
and the average and low Alfven speeds are calculated from OGO-5 observa- 
tions analysed by Chen et al. (1976). As can be seen, the Chiu et el. 
model gives an upper bound to the Alfven speeds. The average OGO-5 
Alfven speeds give the best current Indication of the Alfven speeds 
which are of the order of the beam speed throughout most of the equa- 
torial plasmasphere whose outer boundary is between A and 6 earth radii 
(R e ). Hence V fe < The thermal speeds In Figure 1 are taken from the 

Chiu et al. model. In this case discrepancies between observations and 
the model are unimportant since V^ >> V^ always. 

Another characteristic velocity needing consideration in beam 
modelling is the spread of the beam velocity perpendicular to its direct ion 
of propagation. In the case of the S?S COTV, the spread of the beam 
velocity perpendicular to its direction of propagation is AV^ £ O.Av^ 
(Hanley and Guttman, 1978) as shown in Figure 2 — this corresponds to 
approximately 20° in angular spread and only a 3X loss of effective 
thrust due to the divergence of the beam. This spread is comparable t>< 
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those currently attainable, a.g, a half angle of ^ 15° corresponding to 
only a 2% thrust loss due to divergence Is characteristic of the 30 cm 
Hg systems already developed (D.C. Byers, personal communication). 

Thus, the exhaust of the COTV's may be described as a fast, rapidly 
diverging ion beam. Despite these facts, theories have been developed 
recently (Y. T. Chiu and H. B. Liemohn, private communication) which assume 
that the emerging COTV exhaust will have a very small velocity spread. 

However, the arguments underlying the assumptions of these theories are 
based on laboratory “xperiments whose applicability to the problem of ion 
beam exhaust injection into the plasmasphere remains to be demonstrated. In 
particular, the finite geometries in the laboratory experiments induce 
wall effects which will affect beam spread. In addition if the beam is 
not completely neutralized but contains a substantial net current, 
classical z-pinch (Boyd and Sanderson, 1969) effects will occur which 
will tend to compress the beam. Finally, the laboratory experiments 
were apparently done in parameter regimes with magnetic fields, and 
Alfven speed to beam speed ratios that are not characteristic of the 
plasmasphere. Given these facts, the beam must be regarded as rapidly 
diverging until relevant experiments in space can demonstrate otherwise. 

Due to these beam characteristics, the numerous potential plasma instabilities 
which could take energy from the beam and hence stop it are ineffective. 

This Is due to the fact that the beam and background plasma parameters change 
sufficiently rapidly so as not to allow amplification of instability generated 
wave 8 to significant amplitudes. The interested reader is referred to 
Curtis and Grebowsky (1980a) where a quantitative discussion Is given of the 
criteria needed for instabilities to be a significant beam stopping mechanism. 




Another proposed been stopping theory (Chiu et el., 1979b) models 
the fest Ion been as a slowly moving ion cloud with and 

<< V A< This Is not applicable given the relationship of to 
and V ^ shown in Figure 1. In addition, this ion cloud model assumed 
the beam plasma can be regarded as infinitely conducting. This frozen 
field line concept however is not valid since a realistic model of the 
beam plasma (Ichimaru, 1973) which accounts for both the initial plasma 


turbulence and that generated by the low amplitude plasma wave turbulence 
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carried with the beam gives rapid diffusion times t-X q/D* as shown in 
Figure 3. Note that X^ is the beam ion Debye length and D *, the anomalous 
diffusion coefficient associated with the plasma turbulence, 

D ** 32— (1 + VV W2 ln <» 

where and T g are the ion and electron beam temperatures in energy 
units and m^ and m^, the corresponding masses. The currents resulting 
from the turbulence- Induced anomalous resistivity are insufficient to 
short out the polarization electric field given the range of values of 
AV.. This is true since for a quasi-perpendicular ly propagating ion 


beam the current driven by the polarization electric field, “ — 
is given by 


b 6 
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and since we require AV, > V for the polarisation electric field not 

DC 

to be shorted out we obtain the condition 

Q*B V b (4) 

v c Av b 

for bean propagation. For the expected values of AV^/V^ (Hanley and 
Guttman, 1978) we do not expect the polarisation field to be shorted out 
and hence the beam will propagate. Equation (4) also explains why 
barium releases are often not observed to propagate very far: the initial 
ejection is nearly radial from the explosion point making the spread ot 
the velocity normal to the radial direction small and hence allowing a 
rapid shorting out of the polarization field. In addition, given the 
slower rate of density decrease in the barium releases, plasma instabil- 
ities will be effective in stopping such releases. We emphasize that 
the beam stopping is always a strong function of the initial conditions 
and thus very different behavior may be exhibited by different release 
experiments. Thus, the stopping processes assumed in the ion cloud 
models are not needed to stop either beams or releases. 

Since significant plasma turbulence and an effectively infinite 
plasma conductivity cannot exist simultaneously, ion beam propagation 
models which include both the ion cloud mechanism and plasma instabilities 
(Chiu et ai., 1979b) are not self consistent. Also ion cloud models 
based on the analysis of Scholer (1970) require that the beam's per- 
pendicular velocity spread be aero or very near it, since the model 
treats the beam plasma as an incompressible fluid. Since ion thruster 
exhausts are known to possess large velocity spreads (Hanley and Guttman, 
1978), this is yet another reason why ion cloud models are not applicable 


to the ease of ion bean propagation in tha plasmasphere or in any plasma 
system axe apt those in which tha paraueters change at unrealistically 
slow rates. 

Despite the limitations imposed on bean stopping mechanisms by the 
been velocity characteristics and its finite conductivity, not ell of 
the COTV been pleana will escape the plaanasphare. Since the polarisation 
electric field In posed st the thruster to allow cross-field propagation 
of the been is nonun if or* over the sheath of the bean (see figure 4), 
the plasma In this sheath la lost and deposited on local field lines. 

This model is discussed in detail in Curtis and Grabowsky (1980a). The 
bean sheath formation is a natural consequence of the termination of the 
polarisation field at the sheath's outer edge. A bean sheath loss nodal 
describing this process predicts a deposition of argon ions and hence 
energy in the plasmasphere which is much less than that in models which 
call for ion cloud nechanisms or plasma Instabilities to rapidly stop 
the bean. In Figure 5, a comparison is giver, of the cumulative fractional 
mass loss of an ion bean injected at 1.5 for the ion cloud model and 
the ion beam sheath loss model. The ion cloud procsss stops the beam 
very rapidly whereas all but a few percent of the beau in the ion bean 
sheath loss process escapes. In Figure 6 the integrated difference of 
these two deposition models is shown for the construction of one SPS. 

The ion cloud process results in more than an order of magnitude greater 
energy and number density deposition than that of the sheath loss process. 
The energy spectra of Ar + deposited in the plasmasphere are also dissimilar. 
For the ion cloud model accompanied by a weaker plaaaa instability loss 
process the solid line In Figure '/ gives a qualitative indication of the 
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energy agectr* (Chit: m tl. , 1979b) of the argon lone. In the loa cloud nodal, 
moat of tin energy of Aar + la dlaalpatad in producing Ionospheric curreete ''fftttil 
by the cloud' a fluid line dragging* This process yields the low soar gar peak. 

The higher energy tail and peak just below the injection energy of ^5 IteV would 
be produced by various instability processes. In contrast, the sheath loss 
model shown by the dotted line in Figure 7 results in deposited Ar + with energies 
near the injection energy. 

The ion cloud and ion beam sheath loss models also produce very different 
Ar + plaemasphere densities as a function of time. In the ion cloud model the 
exhaust of a COTV stops rapidly and the time dependent Ar + injection is local. 

In contrast, tha ion ban sheath loss model gives tin dependent Ar + Injection 
that is not localised but rather is extended a number of earth radii away. 

Thus, when the COTV fleet starts transportation from LEO to GEO of SPS construction 
materials, the Ar + COT V exhaust is rapidly deposited throughout the plasmsaphere 
in the ion beam sheath loss model. However, in the ion cloud model the Ar + 
deposition will be localized and move as a front with the COTV flemt. 

To quantify this discussion, we consider the specific case of the Department 
of Energy (1978) scenario for COTV transport of SPS materials using silicon 
photo-voltaic cells. In this scenario, i,_ is projected that there will be IS 
COTV flights per SPS over approximately a 135 day time period. Assuming the 
launches to be equally spaced in time from LEO (i.e., at time intervals of 9 
days) we have done a computer simulation of the time dependent Ar + piasmasphere 
density distribution as a function of time from the first COTV launch from LEO. 

The assumed dominant loss channel for the Ar Is charge exchange with the neutral 
background exosphere. This appears to be the dominant lose channel from our 
previous work (Curtis and Grebowsky, 1980a) and is taken in the computer sim- 
ulation as Che e-folding time for the decay of the injected Ar + . In the calcu- 
lations the relatively small Ar + injection stemming from the COTV fleet return 
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to LEO from GEO is neglected. This is equivalent to assuming that the COTV 
fleet is disassembled at GEO for use in SPS construction rather than being 
reused. The calculation is applicable only to distances within approximately 
four earth radii of LEO because outside of the plasma sphere loss processes other 
than charge exchange may dominate — such as plaama convective loss during substorms. 

In Figure 8 is shown the radial position of a COTV as a function of time 
from LEO departure. As can be seen from the figure, the Injections are greater 
nearer earth (L less than 2.5) than in the more distant plasmasphere. This 
effect is greatly reduced however by the much lower Ar + charge exchange lifetimes 
in the inner region (Curtis and Grebowsky, 1980m). The solid circles indicate 
the radial position of the members of a fleet of COTV's when the first COTV 
reaches GEO assuming a LEO launch approximately <*very 9 days. 

Using the trajectory depicted in Figure 8 the computer simulation -lelds 
the results in Figure 9. The Ar + density distributions are shown at selected 
time intervals of 27, 72 And 133 days for both the ion cloud and ion beam sheath 
loss models. Although the near-LEO Ar + depots cion is greater for the ion cloud 
model than the beam sheath loss model, the sheath model gives a much larger 
spatial extent for the Ar + deposition. Due to the much shorter cltarge exchange 
lifetimes near LEO the total remnant plasmasphere Ar + deposition at day 27 by 
the ion cloud model does >.ot exceed the ion sheath loss deposition by as much 
as the total depoaition difference between the models. Later, et day 72, 
the ion cloud model still is spatially localized near LEO, whereas the 
sheath loss model is characterized by significant depoaition throughout the 
plasmasphere. At day 133, all the materials for the construction of one SFF are 
delivered to GEO and a continuing fleet of COTV's to supply the materials for 
t. _ second SPS are enroute. At this time the A r + densities predicted for the 
ion cloud model are much greater than those for the sheath model. The distinct 
dip in the radial distribution near L • 2.5 is due to the lengthening of lifetimes 
with respect, to charge exchange in the cuter plasmasphere and to the increase in 


L-shell volume In the outer plasmasphere. These effects combine to produce the 
plateau in the Ar + density in the outer plasnasphere. Hie Ar + density plateau 
at 2.5 i L < 5.0 (nominal plasnapauee position) will be the long term signature 
of the As: plaamaspheric deposition. The Inner plasnasphere charge exchange 
lifetimes are sufficiently short so as to give a rapid relaxation of Ar+ back 
to its natural plasnasphere density which is negligible. We note that away from 
LEO the solid angle subtended by the earth decreases fron 2 tt steradians (a full 
hemisphere). Thus at greater distances more of the hot neutral Ar atoms created 
by the charge exchange process will escape into space rather than be recycled 
into the terrestrial atmosphere. 

In Figure 10 is ahov/n a continuation of the computer Simula tier, in tine 
beyond that shorn in Figure 9. The construction materials for a second SPS are 
now being ferried to GEO. We show the resulting time dependent Ar + plasnasphere 
densities as a function of radial distance for tines of 160, 205, and 266 days 
after the first COTV left LEO with the first SPS construction materials. In 
continuing the simulation to the second SPS, it is assumed that the second SPS 
is the last to nave materials transported up to day 266. The plots in Figure 10 
then depict the relaxation of the Ar + densities as the inner source terms are 
removed to successively greater distances — that is, the last member of the COTV 
fleet with construction materials moves outward through the plasmasphere. The 
A r + densities are again plotted for both the ion cloud and the ion beam sheath 
loss models. For both models, the inner plasmasphere densities drop precipi- 
tously after the COTV exhaust source is removed. However, tne outer plasma- 
sphere Ar + density plateau persists due to the long Ar + lifetimes with respect 
to charge exchange in that region. Throughout the plasmasphere during these 
times, the ion cloud model yields much greater Ar + densities than the sheath 
loss model. This has direct implications for environmental effects which scale 
with the magnitude of the Ar + density. 


In summary, from these computer simulations of the Ar + density result- 
ing from COTV exhaust inject Jon into the plasmasphere we have the following 
conclusions : 

1) During the construction of the first SPS, the ion beam sheath loss 
gives an Ar + density enhancement of much greater spatial extent than 
the ion cloud model which produces a perturbation which moves like 

a wave with the COTV through the plasmasphere. The reason for this 
is that in the ion cloud model the beam is stopped in a short dis- 
tance compared to the plasmasphere scale size, whereas in the 
sheath loss model the beam is not stopped. 

2) At later times, after the first COTV reaches GEO, a substantial 
Ar + density plateau develops in the outer plasmasphere. It is 
much greater for the ion cloud model than for the ion beam sheath 
loss model. 

3) The outer plasmasphere Ar + density plateau will tend to persist 
for long time periods given the long charge exchange lifetimes. 

For a simixar reason, the inner plasmasphere Ar + is rapidly depleted. 

In terms of environmental effects, when the results of the computer Simula 

tion arc combined with the differences of the energy spectra of the two models, 

the ion beam sheath loss model may be expected to produce greater possible 

communication perturbations during the Initial transport of materials to GEO 

for the first SPS. This results from the fact that greater numbers of energeti 
+ 

(a. 5 kev) Ar ions will exist over larger spatial extents than in the ion cloud 
model since the beam is not stopped locally. 

Later studies of the Ar + time dependent plasmasphere densities should 
include the consideration of longer term processes involving loss by pitch 
angle scattering and plasmasphere dynamics. This will allow a better 
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assessment of the magnitude of the outer plemsphere Ar + density plateau in a 
Multiple SPS construction scenario. 

In our ion bean sheath loss model «r» have a a sue ad that the Debye length is 
given by an effective tenperature which is proportional to the square of the 
bean spread velocity (Curtis and Grebovsky, 1980a). In essence this implies a 
rapid thotatllsatlon of the perpendicular velocity spread of the bean. Since 
the ion beam emerging as the COTV exhaust is expected to be turbulent upon 
emission due. to the beam creation and charge neutralisation processes, the 
turbulence Implies a rapid randomisation of the perpendicular velocity spread 
component . In particular, in the charge neutralisation process electrons will 
be injected with the ions in the thruster. The differences in velocity in 
directions perpendicular to the ion beam propagation in the electron and ion 
distribution functions may then be expected to drive two-stream type plasma 
instabilities which will generate turbulence and randomise AV^. It should be 
noted that wave amplification process occurs on fixed shells (concentric to 
the beam axis) comoving with the ion beam and hence does not suffer from the 
limitation of wave amplification processes that seek to stop the beam itself 
vuirtis and Grebowsky, 1980a). 

If the thermal isat ion of AV by plasma instabilities is only partially 
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successful, then the effective temperature will be less. Since our ion beam 
sheath loss model Ar + densities scale with the Debye length, they vary as the 
square root of the effective temperature. Partial thermal izat ion could result 
in a factor of ten or more reduction in the Ar density that the model predicts. 

This would make the predicted environmental effects of the ion beam sheath loss 
model even less severe when compared to the massive Ar + deposition envisioned in 
the ion cloud model. To better estimate the degree of beam velocity spread 
thermallsatlon, a better specification of be.im parameters such as emitted icn 
and electron velocity distribution functions is needed than are currently available. 
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The different beam stopping mechanisms can produce very different 
environmental impacts. The sheath loss model predicts a large injection 
of energetic anisotropic argon ions which will drive plasma instabilities 
which may produce sufficient scintillation to Impair radio communications 
with geosynchronous satellites. The partial depletion by precipitation 
of the energetic ion belts surrounding the earth is also possible due to 
the pitch angle scattering caused by Ar + turbulence. A detailed description 
of potential environmental effects using the sheath loss model is given 
in Curtis and Grebowsky (1980b). Cold argon ions (T ^ 1 eV) are produced 
in the sheath loss model only via the loss of energy by plasma instability 
mechanisms and electron coulomb scattering. Since during the energy 
degradation processes, Ar + will be lost by charge exchange and precipi- 
tation, the amount of cold Ar + plasma from the sheath loss mechanism 
will be much less than from the ion cloud mechanism. The environmental 
effects due to cold Ar + would be greatly reduced in the sheath loss 
picture as well as those effects due to ionospheric currents (Chiu et 
al. , 1979b). Specifically, with the greatly reduced numbers of cold Ar + 
expected in the ion sheath model, the dosage enhancement of trapped 
relativistic electrons predicted by Chiu et al. (1979b) would not occur. 

The extension of the radiation belt environment out to the geosynchronous 
region envisioned by Cladis and Davidson (1980) using the Chiu et al. 

(1979b) model would also not occur since it requires large cold ion 
deposition. The artificial ionospheric current predicted by the Chiu et 
al. model relies exclusively on the ion cloud model's frozen- in field 
lines being dragged. Without the ^ery doubtful "frozen-in picture", the 
Chiu et al predicted powerline tripping and pineline corrosion will not 
occur. To claim that results of the ion cloud and beam sheath loss 
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medals fare eke sue because the final state In both cases is to leave 
cold Ar + in the plasma sphere Is meaningless. The radically different 
energy deposition spectra shown in figure 7 yield very different plasma- 
sphere modifications before reaching their final state (cold Ar + ) by 
various energy degradation processes. 

Relevant Experiments 

To establish a convincing final assessment of the beam propagation 
physics, experimental tests are required to overcome the credibility 
problems. As Roger Bacon said in his work Opus Tertium. "The strongest 
of arguments prove nothing so long as the conclusions are not verified 
by experience." This is true In plasma physics in particular where 
experiment characteristically leads theory. From this viewpoint, a very 
good experimental test would be a Space Si.uttle-born ion beam experi- 
ment. Remote in situ measurements of the ion beam characteristics could 
be done by a free flying satellite launched prior to the experiment by the 
space shuttle. This could be a scaled down COTV Ar + thruster with power levels 
of about a kilowatt and a nozzle diameter of a few centimeters rather 
than the COTV values of a megawatt and a meter respectively. The development 
of suitable advances in ion thrusters for such an experiment does not sets 
to pose major technical problems but is rather limited by the willingness 
to expend the necessary funding (Hudson and Finke, 19d0) . The other 
beam parameters could be the same as for a COTV. The down sized COTV 
thrusters would then have the same exhaust beam speed, perpendicular velocity 
spread, density and ion mass as the projected full sized COTV thrusters. The 
required power levels could be within the limits of the planned 20kW orbiter 
integral solar array or the 6kW orbiter mounted array (Franklin, 1980). 

We emphasize that experimental conditions which do not closely 
correspond to those of the COTV icn beam exhaust will not give a meaning- 
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ful assessment of the exhaust physics. In particular, In order to be 
applicable, the Ion mass, the Initial beam density, the beam velocity 
and Its perpendicular velocity spread as well as, the background plasma 
density, composition and magnetic field used in experiments must be near 
those expected in the case of COTV Ion beam Injection. Only in space 
can these conditions be well satisfied. 


Mitigation Strategies 

If the perpendicular velocity spread of the ion beam exhaust can 
actually be reduced, despite current engineering problems, to where 


A V b 


«1, then plasma instabilities will become important in stopping 


the beam as discussed in Curtis and Grebcwsky, (1980a). However this 

e 

effect can be reduced by increasing AV^/y and thus increasing n^ n 

b b 

and not allowing large amplification of instabilities. Increasing 


AV, 


b/V. 


may also be necessary to stop the turbulence driven currents from 


shorting out and thus stopping the beam. This will, however, reduce 
the overall efficiency of the propulsion system since AV^ does not 
contribute in lifting cargo to higher orbits. When argon was chosen as 
the propellant the choice was based on its availability and the comblnal U 
of high specific impulse and high thrust that using an intermediate mass 
ion allowed. Use of either a lighter or heavier mass ion would sacrifice 
these overall efficiencies. For a given beam energy, the beam velocity 
scales as m^ for a specific ion mass. Instability processes involving 
the beam are characteristically proportional to the beam ion plasma 
frequency u (A” n e an< * t * ie * on gy rc_ f re <l uenc y ^ “ /m^c. 

-1 Ja 

Thus n,* m. and oj « m, . Plasma instabilities for a fixed beam 
i i pi i 

energy can thus be suppressed by choosing a heavier mass ion since the 
allowed length for instability amplification is proportional to 1 

where the growth rate <d^ a or ui^a Choice of a heavier ion will 
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however dtcraue the beam speed for a fixed beam 
also be necessary to increase in order to suppress instabilities 
after the rapid beam expansion phase. In the later phase of beam 
propagation the restriction on instability growth lengths are given by 
the length scalea of the convective derivatives discussed in Curtis and 
Grebovsky (1980a) that arc related to inhomogenietles in the plasmasphere 
and are proportional to V^. Increasing may be necessary even without 
going to a heavier ion if plasma instabilities are effective in stopping 
the beam in the post-rapid expansion phase of its propagation. 

Conclusions 

Understanding the correct physics of propagation of the ion beam 
emitted by the CQTV's envisioned for SPS is critical to assessing what 
environmental effects may occur. In the two models discussed here the 
underlying physics is very different as are the predicted environmental 
effects. We believe that our theory, which models the ion beam's propa- 
gation by a sheath loss model, is much more in accord with the present 
understanding of plasma physics than is the ion cloud model which assumes 
concepts such as "frozen- in" field lines without consideration of the 
restricted range in which the assumptions are valid. Our theory is also 
self consistent - in contrast, for the reasons mentioned above, the ion 
cloud model is not and (as applied to the COTV ion beam problem) often 
calls for the simultaneous operation of mutually exclusive processes. 

The major environmental effect that we foresee from our modelling of the 
COTV ion beams is communication impairment du^ to scintillations caused 
by ionospheric and plasmasphorlc tensity irregularities . These irregu- 
larities are produced by the hot anisotropic argon ions deposited by the 
beam. The more spectacular effects predicted by the ion cloud model 


energy E^, so it may 
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will not occur or will be reduced in magnitude by many ordera of magnitude 
due to the reduced deposition of beam ions deduced in the beam sheath 
model. Should technological advances require modifications to the system 
design of the COTV thrusters, the beau model developed by us can be used 
to determine how best to minimize the environmental impact through judicious 
selection of the free beam parameters. 

We strongly feel that that appropriate space, as opposed to laboratory, 
experiments are needed. Since the capability to perform these studies will 
not be obtained until about the mid- 80 ' s, we recommend that additional 
theoretical work be done in order to lay a firm foundation for the questions 
to be addressed by the experiments. Specific theoretical work suggested 
includes both a more refined modelling of the ion beam propagation in *"he 
plasmasphere, of the deposition of beam ions in the plasmasphere, and of 
the time dependent overall plasmasphere response to this input. 
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FIGURE CAPTIONS 


Figure 1 . Comparison of the characteristic Alfven speeds and thermal ion speeds 
with the projected COTV argon ion exhaust speed as a function of dis- 
tance in earth radii in the equatorial plane. 


Figure 2. Diagram of expanding argon ion beam emitted by a COTV. The width of 
the beam sheath, the Debye length Ap, is shown in an enlargement . 

Figure 3. Diffusion time of magnetic field through the turbulent ion beam 

plasma. The field lines will rapidly slip through the beam rather 
than being "frozen-in" . 


Figure 4. Gradient of beam speed across the beam sheath due to the drop off of 
the pola-'zation electric field ove. the Debye length X_. This gives 
a continual shedding of the beam sheath and hence yields the deposition 
of energetic, anisotropic argon ions in the plasmasphere. 


Figure 5. Comparison of fraction of total mass of an ion beam deposited in the 
plasmasphere as a function of radial distance when the COTV is at 
1.5 earth radii. Two beam models are shown: ion cloud and ion beam 
sheath loss. 

Figure 6. Comparison of total number density of Ar + deposited in the plasmasphere 
for each SPS built. The ion cloud and ion beam sheath loss models are 
compared . 

Figure 7. Comparison of the energy spectra of the Ar + deposited in the plasma- 
sphere by the ion cloud and ion beam sheath loss models. 

Figure 8. Radial position of a COTV carrying SPS construction materials 
as a function of time. The points on the curve indicate the 
radial positions of the members of a COTV fleet if launching 
occurs equally spread in time. 

Figure 9. Time dependent deposition of Ar + in the plasmasphere for a 
COTV fleet carrying SPS construction materials from LEO to 
GEO. The silicon scenario SPS parameters are used here 
(Department of Energy, 1978) and the COTV's are assumed to 
depart from LEO in an equally spaced time sequence. Deposition 
for ion cloud and ion beam sheath models — the solid and dashed 
curves respectively — are compared. 

Figure 10 Continuation of Figure 9 for the COTV transport of construction 
materials for a second SPS. The launchings from LEO are again 
assumed to be equally spaced in time and start with the first 
time interval after the last COTV for the first SPS is launched. 
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BEAM DEBYE SHEATH LOSS 


O 

O II 

•• CD 
III > 

h- 

z 

ui 

5 

LU 

CD 

QC 

z 

LU 




” I 
LU H 

iSg 

oc □ 

0 

lu w 

Q >: 
CD 

RS 

DO 

cog 

o< 

_J 

1 
H 

IS 

S 


CL 

LU 


LU 



FIGURE 4 


t 


(3) SHEATH LOST EVERY GYROPERIOD 


CUMULATIVE FRACTIONAL LOSS 


DENSITY LOSS OF ION BEAM 



FIGURE 5 



NUMBER DENSITY (NO/CC) 





GEOCENTRIC DISTANCE (R e ) 


FIGURE 6 



NORMALIZED DENSITY PER SPS 


COMPARISON OF ION CLOUD 
AND ION BEAM SHEATH LOSS 
MODEL DEPOSITED ENERGY SPECTRA 


ION CLOUD 


10-1 


10-2 


ION BEAM 
SHEATH LOSS 


IO - 3 

10 - 


ENERGY (EV) 


FIGURE 7 





i*. m un « m n 

i 3 n/iin 

FIGURE 8 












